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Abstract 

Phase transitions in hot and dense matter and the in-medium behavior of pseudoscalar mesons 

T]') are investigated, in the framework of the three flavor Nambu- 
Jona-Lasinio model, including the 't Hooft interaction, which breaks the Ua{^) symmetry. Three 
different scenarios are considered: zero density and finite temperature, zero temperature and finite 
density in quark matter with different degrees of strangeness, and finite temperature and density. 
At T = 0, the role of strange valence quarks in the medium is discussed, in connection with the 
phase transition and the mesonic behavior. It is found that the appearance of strange quarks, 
above certain densities, leads to meaningful changes in different observables, especially in matter 
with (3 -equilibrium. The behavior of mesons in the T — p plane is analyzed in connection with 
possible signatures of restoration of symmetries. 
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I. INTRODUCTION 



Understanding the behavior of matter under extreme conditions is nowadays a challenge 
in the physics of strong interactions. Different regions of the QCD phase diagram are object 
of interest and major theoretical and experimental efforts have been dedicated to the physics 
of relativistic heavy-ion col·lisions, looking for signatures of the quark gluon plasma [QGP] 
. Special attention has also been given to neutron stars which are a natural laboratory 
to study matter at high-densities. 

There are indications from Lattice QCD that the transition from the hadronic phase to 
the quark gluon plasma is probably associated to the transition from the Nambu-Goldstone 
realization of chiral symmetry to the Wigner-Weyl phase. While the phase transition with 
finite Chemical potential and zero temperature is expected to be first-order, at zero chemical 
potential and finite temperature there will be a smooth crossover. Experimental and theo- 
retical efforts have been done in order to explore the ^ — T phase boundary. Recent Lattice 
results indicate a critical "endpoint", connecting the first-order phase transition with the 
crossover region, at Te = 160 ± 35MeV, /i^; = 725 ± 35MeV 4]. Understanding the results 
of experiments at BNL Q and CERN P| provides a natural motivation for these studies. 

The rich content of the QCD phase diagram has been recently explored in the direction of 
high-density and cold matter, that can exist in neutron stars, where a " color-fiavor locking" 
[CFL] phase, exhibiting a variety of interesting physics is expected to occur 

Restoration of symmetries and deconfinement are expected to occur at high-density and 
/or temperature. In this concern, the study of observables of pseudoscalar mesons is particu- 
larly important. Since the origin of these mesons is associated to phenomena of spontaneous 
and explícit symmetry breaking, its in-medium behavior is expected to carry relevant signs 
for possible restoration of symmetries. On the classical level and in the chiral limit, the 
QCD Lagrangian has two chiral symmetries, the SU{Nf) and the ^7^(1). This would imply 
the existence of nine Goldstone bosons, for Nf = 3, and, in order to have mesons with 
finite mass, chiral symmetry must be explicitly broken ab initio by giving current masses 
to the quarks. However, in nature there are only eight light pseudoscalar mesons, the octet 
(vr, K, T]); the rj' has a mass too large to be considered a remnant of a Goldstone boson, so 
its mass must have a different origin. In fact, the Ua{Í) symmetry is not a real symmetry of 
QCD, since it is broken at the quantum level, as pointed out by Weinberg The breaking 
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of the í/a(1) symmetry can be described semiclassically by instantons, which has the effect 
of giving a mass to Vj of about 1 GeV. On the other side, this í/4(l) anomaly causes flavor 
mixing, which has the effect of hfting the degeneracy between r\ and tt^. So a percentage of 
the ?7 mass comes from the U symmetry breaking. 

So far as the restoration of symmetries is concerned, there are two possible scenarios 
only ^[/(S) chiral symmetry is restored or both, S'f/(3) and ?7a(1) symmetries, are restored. 
The behavior of r/' in medium or of related observables, like the topological susceptibility 
P], might help to decide between these scenarios. A decrease of the rj mass could lead to 
the increase of the r( production cross section, as compared to that for coUisionsJ ^ . 

Strange quark matter [SQM] has attracted a lot of interest since the suggestion IjJ that 
it could be the absolute ground state of matter. Stable SQM in /3 -equilibrium is expected 
to exist in the interior of neutron stars or even be the constituent matter of highly bound 
compact stars ("strange quark stars"). Lumps of SQM, the strangelets, might also be formed 
in earlier stage of heavy-ion col·lisions (in this case /3 -equilibrium may not be achieved). 
Experiments of ultrarelativistic heavy-ion col·lisions at BNL and CERN are proposed to 



search for strangelets 



but up to now there is no evidence of such objects 



Several studies on the behavior of matter with different strangeness fractions are known (see 



and references therein). 
An interesting pro blem in flavor asymmetric matter is the behavior of the charge mul- 
tiplets of mesons [lllllSl· These charge multiplets, that are degenerated in vacuum or in 
symmetric matter {p^ = Pd = Ps), are expected to have a splitting in flavor asymmetric 
matter. In particular, the masses of kaons (antikaons) would increase (decrease) with den- 
sity. A similar effect would occur for tt^ and tt~^ in neutron matter. A slight raising of the 



mass and a lowering of the K mass 
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2G 
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23| seems to be compatible with 



24| and with the results of KaoS and FOPI coUabora- 



the analysis of data on kaonic àtoms 

n 

tions at GSI [2^. The driving mechanism for the mass splitting is attributed mainly to the 
selective effects of the Pauli principle, although, in the case of K~ , the interaction with the 
A(1405) resonance plays a significant role as well, in the low-density regime 0|. 

Effective quark models are useful tools to explore the behavior of matter at high-densities 
or temperatures. Nambu-Jona-Lasinio [NJL] [26] type models have been extensively used 
over the past years to describe low-energy features of hadrons and also to investigate restora- 
tion of chiral symmetry with temperature or density 
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34, 



35 



36|. 



The NJL model is an effective quark model where the gluonic degrees of freedom are 
supposed to be integrated out and, besides its simplicity, has the advantage of incorporating 
important symmetries of QCD, namely chiral symmetry. Since the model has no confining 
mechanism, several drawbacks are well known. It should be noticed, for instance, that the 
rj mass is very close to the qq threshold and, depending on the parameterization, it can 
be above this threshold. In this case, this meson is described, even at zero temperature 
and density, as qq resonance, which would have the unphysical decay in qq pairs, and the 
definition of its mass is unsatisfactory. 

The behavior of S't/(3) pseudoscalar mesons in hot matter has been studied within the 

34, [3^. Different studies have been devoted to the 

or asymmetric 



framework of NJL model in 



behavior of pions and kaons at finite density in flavor symmetric 



matter 
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27, 



A model aiming at describing hadronic behavior in the medium should account for the 
great variety of particle-hole excitations that the medium can exhibit, some of them with the 
same quantum numbers as the hadrons under study 0, l^l • Particle-hole excitations with 
the same quantum numbers of kaons have been discussed in jïsl. It has been shown, within 



the framework of NJL models, that low-energy pseudosca^ 



ar modes, which are excitations of 
29|. Such studies were carried out 
0), for charged kaons, and neutron 



the Fermi sea, occur in flavor asymmetric media 27, 
in quark matter simulating nuclear matter i^p^ = pd , ps 
matter without (3 -equilibrium, for charged pions. The role of the 't Hooft interaction was 
not taken into account in the last case. The combined effect of density and temperature, as 



19,^ 



well as the effect of vector interaction, was discussed for the case of charged kaons 
Only densities below ~ 3po were considered. The high-density region (p„ > p'^ = 2.25po) of 
quark matter simulating neutron matter in weak equilibrium, having in mind the study of the 
behavior of kaonic (charged and neutral) and pionic (charged) excitations was investigated 
in hTf, and the behavior of neutral pseudoscalar mesons in hot and dense matter was 

n 

investigated in j38| . 

This paper is devoted to investigating the phase transition in hot and dense matter, 
and the in-medium behavior of the pseudoscalar mesons, in the framework of the SU{3) 
Nambu-Jona-Lasinio model, including the 't Hooft interaction, having in mind to look for 
manifest ations of restoration of symmetries and to discuss the role of the strangeness degree 
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of freedom. 

We present the model and formalism in the vacuum (Section|n|) and at finite density and 
temperature ÍSection We will present and discuss our results for the phase transition 
and meson behavior in three scenarios: finite temperature and zero density ÍSection II V|) : 
zero temperature and finite density in quark matter with different degrees of strangeness, 
with and without f3 -equilibrium (Section|V|). Finally, the meson properties in hot and dense 
matter are investigated in Section IVll In Section [VIII we draw our conclusions. 



II. MODEL AND FORMALISM 



The NJL model with three quarks can be described by the Lagrangian 



a=0 



+ go 



det[g(l + 75)g] + det[g(l - 75)5] 



(1) 



Here q = {u, d, s) is the quark field with three fiavors, Nf = 3, and three colors, Nc = 3. 
A'^ are the Gell-Mann matrices, a = 0, 1, . . . , 8, A° = -y/ï I. The explícit symmetry breaking 
part in contains the current quark masses rh = diag(m„, m^^, m^). The last term in is 
the lowest six-quark dimensional operator and it has the SUlÍ'Í) ® SUr{3) invariance but 
breaks the Ua_{Í) symmetry. This term is a refiection of the axial anomaly in QCD. For 
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general reviews on the three fiavor version of the NJL model see 

FoUowing a Standard hadronization procedure, the foUowing effective action is obtained: 

1 



Weff[ip,a] 
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—íTt In 



-m + aaX" + (Í75)(<^aA") 



(2) 



The notation Tr stands for the trace operation over discrete indices {Nf and Nc) and inte- 
gration over momentum. The fields a°' and y?" are the scalar and pseudoscalar meson nonets 
and Sab , Pab are projectors defined in the Appendix A. 

The first variation of the action Q leads to the gap equations, 



Mi = mi- 2gg < qiqi > -2g^ < qjqj >< qkqt > 



(3) 



with k = u,d,s cyclic. Mj are the constituent quark masses and the quark condensates 
are given by: < (jiqi >= —iTr[Si{p)] , Si{p) being the quark Green function. 



To calculate the meson mass spectrum, we expand the effective action (j21) over meson 
fields. Keeping the pseudoscalar mesons only, we have the effective meson action 



(2) 



--09" 

2^ 



p-,' - n,,(P) 



(4) 



with I{ab{P) being the polarization operator (see Appendix A). The expression in square 
brackets in (ji} is the inverse nonnormahzed meson propagator D^^{P). The pseudoscalar 
meson masses are obtained from the condition (1 — PijW^{PQ = M,P = 0)) = 0. For the 
nondiagonal mesons tt , K, the polarization operator takes the form: 



U^Po) = 4 ((/i + n) - [P^ - (M, - M,)^] I^Po 



(5) 



where the integrals are given in the Appendix A. 

The quark-meson coupling constants are evaluated as 

1 d 

- IV-'iPn) 

\P0=M 



9m; 



(6) 



2M dPo 

where M is the mass of the bound state containing quark fiavors 

To consider the diagonal mesons vr", rj and rj' we take into account the matrix structure 
of the propagator in (0)). In the basis oí tt^ — rj — t]' system we write the projector Pab and 
the polarization operator Hat as matrices: 



ab 



^ P33 P30 P38 ^ 

-P03 -Poo -P08 



\ P83 PsO Psi 



( Hss nso Has ^ 



and n, 



ah 



I 



Hos 

Hss IIso Hgs j 

The nondiagonal matrix elements P30 = "^5'd(< Qu^u > — < qdQd >), P 



(7) 



38 



733 D 



< 



ququ > - < QdQd >), ^30 = y^2/3[J„„(Po) - Jdd{Po)] and Uss = l/V5[Juu{Po) - Jdd{Po)] 
correspond ton^ — r] and n^ — rj' mixing. In the case < Qulu >=< qdld >, the 7r° is decoupled 
from the r] — rj' and these elements vanish. The specific form of the nonvanishing elements 
of those matrices may be found in the Appendix A. 
Defining the orthogonal matrix 

0=1 V (8) 



COS0 — sin^ 
smO cosO 
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we may find the r] — i]' mixing angle 9 via the condition to diagonalize D^^{P) as 
0-^D^^{P)0 = diag(D-^(P),L'^/(P)). To find the masses and M^/, we use the in- 
verse propagators 



D-\P) = {A + C)- V(C - Ay + 4i?2 (9) 



D~,\P) = {A + C) + ^iC-Ay + 4B^ (10) 

with A = P88-Anoo(P),C = Poo-AUssiP),B = -{Pos + Aïl^siP)) and A = PqoAs-í'os· 
In the rest frame, D-\Po = M^, P = 0) = 0, D-,\Po = M^/, P = 0) = 0. The mixing angle 
9 can be calculated by tan26' = 2B/{A — C). The coupling constants are determined by (jü)). 
Note that from (jü)) we have the quantities gQq.gsq and gori',g8r]', from which we can obtain 



gr){·q')üu and grf{·q')ss, using Standard expressions (for more detalls see for example j3J,|39|). 

The model is fixed by the coupling constants gs,gD in the Lagrangian ((T)), the cutoff 
parameter A which regularizes momentum space integrals !{ and l2{P) and the current 
quark masses mj. For our numerical calculations we use the parameter set js^ : 
mu = md = 5.5 MeV, m, = 140.7 MeV, gsA^ = 3.67, goA^ = -12.36 and A = 602.3 MeV, 
that has been determined by fixing the vàlues 

= 135.0 MeV, Mk = 497.7 MeV, = 92.4 MeV and M^, = 960.8 MeV. We also have 

= 514.8 MeV, ^(M^) = -5.8°, g^üu = 2.40, g^^ss = -3.91, 
9{M^,) = —43.6°, grj'üu = 2.69, grj'ss = —0.54. For the quark condensates we have: 
< q^q^ >=< q^q^ >= -(241.9 MeV)^ and < >= -(257.7 MeV)^ and Mu = Ma = 
367.7 MeV, Mg = 549.5 MeV, for the constituent quark masses. 

III. THE MODEL AT FINITE TEMPERATURE AND CHEMICAL POTENTIAL 

Now we generalize the NJL model to the finite temperature and chemical potential case. 



It can be done by the substitution (see |40|) 



d'p l f d^P ^ 



(27r)4 -zf3J (27r)3^ 

where /3 = 1/T, T is the temperature, and the sum is done over Matsubara freqüències 
üün = (2rï, + l)7rT, n = 0, ±1, ±2, . . ., so that po — > iuJn + with a chemical potential /i. 
Instead of integration over pq, we have now the sum over Matsubara freqüències which can 
be evaluated as 
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f{zm)) Res[h{uj 



where f{z) and f{z) are the Fermi distribution functions for quarks and antiquarks, 



(12) 



As 1 — /(z) = f{—z), we introduce, for convenience, the Fermi distribution functions of the 
positive (negative) energy state of the ith quark: 

The integrals !{ , I2 (-P) that enter in the expressions of the propagators depends now 
on the temperature T and on the chemical potentials, in a Standard way (see Appendix 
A). Having these integrals, we can investigate the phase transition and meson properties in 
hot and dense matter. First, we analyze the mesonic behavior at finite temperature and 



vanishing chemical potentials. Although such a study was already performed in [31|,|3J], it 
is pertinent to present the results here for the sake of comparison with our new findings at 
finite temperature and density, to be presented in next sections. 

Fig. 1 shows the temperature dependence of vr, K, rj and 1]' meson masses, as well as 
of 2Mjj and M„ + Mg at /i = 0. One can see that, at low temperature, the masses of 
mesons (except the ï^'-meson that is always unbound) are lower than the masses of their 
constituents. In this case the integrals I2 are real. The crossing of the vr and 77 lines with 
2Mu and the K line with M„ + Mg indicates the respective Mott transition temperature 
Tm for these mesons. Mott transition comes from the fact that mesons are not elementary 
objects but are composed states of quarks, and is defined by the transition from a bound 
state to a resonance in the continuum of unbound states. Above the Mott temperature we 
have taken into account the imaginary parts of the integrals I2 and used a finite width 



approximation [3^|35[. For our set of parameters, one can see that Mott temperatures for 
T] and vr mesons are: Tm-tj = 180 MeV and Tmtt = 212 MeV. The vr and K mesons become 
unbound at approximately the same temperature. 
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FIG. 1: Temperature dependence of the pion (sòlid line), kaon (dotted line), rj (dashed line) and rj' 
(short-dashed line) masses. The curves 2Mu and M^ + Mg (dot-dashed line) show the temperature 
dependence of the quark thresholds. The respective Mott temperatures are: T/v/tt — Tmk = 212 
MeV and Tm^ = 180 MeV. 

IV. THE PHASE TRANSITION 

The nature of the phase transition in NJL type models at finite T and/ or /i has been 



discussed by different authors 
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4l| . At zero density and finite temperature 



there is a smooth crossover, at nonzero densities different situations can occur. We will 
discuss this problem by analyzing the behavior of the pressure and of the energy per particle 
as functions of the baryonic density. 

The equilibrium state may be determined as the point where the thermodynamical po- 
tential takes the minimum with the quark condensates < g^gj > as variational parameters. 
The baryonic thermodynamic potential has the foUowing form: 

n{p,T)=E-TS- l·^^Ni^ (15) 

i=u,d,s 

where E, Sa.ndNi are, respectively, the internal energy, the entropy and the number of 
partides of the ith quark, that are given by the foUowing expressions: 

E = -—^V Y \j V dv ^ {n^ - O ^(A - p ) ^ 

i=u,d,s L * J 



-9s ({qiqi))^ -'^gD{uu){dd){ss) , 

i=u,d,s 



(16) 



S - -^V E / P^rfp ^(A^ - P^) 

i=u,d,s 

X I [nt In n+ + (1 - n+) ln(l - n+)] + [n+ ^ nr] } , 



(17) 



= / P'^P í^^- + - l) ^(^' - P') ■ (18) 
V is the volume of system and the quark density is determined by the relation pi = Ni/V. 

We define the pressure and the energy density such that their vàlues are zero in the 
vacuum: 

p(p,r) = -i[Q(p,r)-Q(o,r)], (i9) 

e{p,T)^^[E{p,T)-E{0,T)]. (20) 
A. Chiral phase transition at zero temperature 

We start by analyzing the behavior of quark matter at zero temperature and, in order 
to discuss the role of the strangeness degree of freedom, we will consider matter with and 
without (3 -equihbrium and, in the last case, we assume different fractions of strange quarks. 

In "neutron" matter in chemical equihbrium, maintained by weak interactions and with 
charge neutrahty, the foUowing constraints on the chemical potentials and densities of quarks 
and electrons should be imposed: 

Aid = Ais = Aí« + Aíe, (21) 



■^Pu- :^{pd + Ps) - Pe = 0, (22) 

with 

p, = ^(/.^ - Mlf'^d{i4 - Ml) andpe = A^^/Stt^. (23) 

We should note that if the chemical potential of the electron exceeds the rest mass of the 
muon (//g > Mjj) it becomes energetically favorable for an electron to decay into a pL~ via the 
weak process e~ — > + ï?^ + and we can have a Fermi sea of degenerate negative muons. 
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Consequently, the condition for charge neutrality would be fpw — |(pd + Ps) — Pe — = 0. 
However we have found p^"^ = 95.7 MeV < M^. So, we can neglect the muon contribution. 

The thermodynamic potential, pressure and energy defined in (fTÜ|) and piH) have to 
be modified in order to include the contribution of electrons. This leads to the following 
expressions for the pressure and energy density (see 

P = -e(p,0) + J2 /^«P^ (24) 

i=u,d,s 

^ = e(p,0)+^. (25) 
Therefore the zeros of the pressure give the following expression for the energy density: 

^ = J2 l^iPi + l^epe- (26) 

i=u,d,s 

Defining the baryonic matter density as p„ = ^{pu + Pd + Ps) and using the conditions (^1]) 
and ()22|) it is straightforward to show that the energy per baryon at the zeros of the pressure 
takes the form 

— = p„ + 2pd. (27) 

Pn 

The pressure has a zero at p„ = 0, the energy per baryon at that point being Mu + 2M(i (since 
in vacuum the masses of the quarks m, d are equal we will from now on denote this quantity 
as 3M^). If there is another zero of the pressure, at p„ 7^ 0, that corresponds to a minimum 
of the energy, the criterion for stability of the system at that point is pu + 2pd < 3M„. Let's 
now analyze our results for the pressure and energy per baryon, that are plotted in Fig. 2 
a)-b) as functions of p„/po, where po = 0.17 fm~^ is the normal nuclear matter density. The 
pressure has three zeros, respectively at p„ = , 0.45po , 2.25po, that correspond to extrema 
of the energy per particle. For p„ < 0.2po the pressure and compressibility are positive, so 
the system could exist in uniform gas phase, but it will not survive indefinitely, since the 
zero density state is energetically favored; for 0.2po < p„ < 0.45po the system is unstable 
since the compressibility is negative, in fact p„ = 0.45po corresponds to a maximum of the 
energy per particle; for 0.45po < p„ < 2.25po, the pressure is negative, and the third zero 
of the pressure, p = 2.25po, corresponds to an absolute minimum of the energy. In fact, at 
that point p„ + 2pd = 1099.4 MeV, lower then 3M^ = 1102.9 MeV. Above p„ = 2.25po, that 
we define as pcr, we have again a uniform gas phase. The phase transition described in this 
model is a first-order one since there is a region where the pressure and/or compressibility are 
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FIG. 2: Energy per baryon number (a)) and pressure (b)) as a function of density. Sòlid line: 
r = 0, dashed line: T = 56 MeV. 

negative and, given that at the critical density /i„+2yUrf < 3M^ (we notice that the satisfaction 
of this condition is dependent on the parameter choice), the behavior described allows the 
interpretation that the uniform gas phase will break up into stable droplets of nonstrange 



quarks with partially restored chiral s 



nontrivial vacuum ísee also 



m 
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ynimetry and density pjf = 2.25po5 surrounded by a 

m). 



We will now discuss cases of matter without /3-equilibrium. As we will show, the results 
are qualitatively similar to those discussed above, although some specific aspects deserve a 
closer analysis. We consider three cases: Case I - "neutron matter" without strangeness, 
{pd = 2pu , Ps = 0); Case II - matter with equal chemical potentials {p = Pd = l^u = /^s) 
with isospin symmetry, pu = Pd, Ps = ~2 (/U^ — M'^Y^'^6{p'^ — Mf); Case III - matter entirely 
fiavor symmetric {pd = Pu = Ps)- 
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Let us make some comments concerning Case III, which intuitively seems the less natural 
scenario. AlthouErh rather schematic, this case simulates a situation where we can explore 
the Kypothes. of ab.o.utelv stable SQKl It Kas been .,„ed Q Q that SQM ,„a. „„ly 
be stable if it has a large fraction of strange quarks {ps ^ Pu ^ Pd)- The speculations on 
the stability of SQM are supported by the following observations: the weak decay of strange 
quarks into nonstrange quarks can be suppressed or even forbidden due to Pauli blocking, 
and, in addition, the inclusion of a new flavor degree of freedom allows for a larger decrease 
of strange quark mass which can produce a sizable binding energy. As we will see, Case III 
confirms this tendency when compared with, for instance, Case I. 

So, in spite of the simplicity of our assumptions, we think that the analysis of these types 
of matter could be useful as a guideline to understand the role of the strangeness degree of 
freedom in matter that is supposed to exist in neutron stars and in matter that might be 
formed in an early stage of heavy-ion col·lisions. This will also be relevant for the discussion 
of the mesonic behavior that will be presented in next section. 

We notice that in all cases there is an absolute minimum of the energy per particle at 
nonzero density and zero pressure, lower than the vacuum constituent quark masses, so 
we have a first-order phase transition with the formation of quark droplets (see Fig. 3). 
However, this minimum of the energy is always higher than 930 MeV, which is the energy 
per baryon number in atòmic nuclei, and therefore we do not find absolutely stable quark 
matter. The lower energy per baryon is for Case II, and the higher binding energy, compared 
to the vacuum constituent quark masses, is for Case III. The minimum of the energy occurs 
at pn = 2.25po in matter in (3 -equilibrium and in Case I, at p„ = 2.33po for Case II, and 
at Pn = 3.50po for Case III. Only in the last case we find stable SQM, since in the other 
cases the minimum occurs in the region where strange quarks do not exist (see Fig. 4, right 
panel). 

Finally we add some considerations about the restoration of chiral symmetry. When 
chiral symmetry is broken ah initio, one can talk about restoration of chiral symmetry 
if the constituent quark masses drop to the current quark masses, indicating a transition 
from the chirally broken to approximately chirally symmetric phase. This would happen at 
asymptotic densities, when the Fermi momentum equals the cutoff. However, at densities 
considerably lower, one can see, or not, a clear tendency to the restoration of chiral symmetry 
depending on the quark sector considered and on the strangeness fraction ^ . To illustrate 
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FIG. 3: Energy per baryon number for all types of quark matter considered at T = 0. 

this point we plot in Fig. 4 the dynamical quark masses and chemical potentials as functions 
of pnl po- The nonstrange quark masses decrease sharply in all cases, reflecting the quick 
restoration of chiral symmetry in this sector and, as cxpcctcd, the behavior of the strange 
quark mass depends strongly on the amount of strange quarks in the medium. In all cases, 
except Case III, the mass decrease up to ~ 2po is due to the 't Hooft contribution in the gap 
equations. In Case I it remains constant afterwords, since strangeness is put equal to zero 
by hand; in Case II and in (5 -equilibrium, at densities above 5po and 3.8po, respectively, the 
mass becomes lower than the chemical potential and a more pronounced decrease is observed. 
These densities are the onscts for the appearancc of strange valence quarks, which become 
more important as the density increases. In Case III, there are always strange valence 
quarks present, so the strange quark mass decreases more strongly, although, even in this 
case, is still away from the strange current quark mass. We observe that in matter in 
(5 -equilibrium there is a tendency to the restoration of flavor symmetry, as the density 
increases: the chemical potentials — /ig and /x^ become closer and /ig decreases. 
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FIG. 4: Constituent quark masses Mu, M^, Mg and chemical potentials fiu, IJ-d, fJ-e as functions of 
density at T = 0. Left panel: without /3 -equilibrium {Mu = in Cases II and III). Right panel: 
with (3 -equilibrium. 

B. Chiral phase transition at finite temperature and density 

Now we discuss the phase transition in hot and dense matter and we consider only matter 
in (3 -equilibrium (the other cases are qualitatively similar). Since for very low temperatures 
the absolute minimum of the energy turns to be at zero density, the phase transition is still 
íirst-order but the system is unstable against expansion. With increasing temperature, we 
will have a crossover above T > T'c/ = 56 MeV. The point T = 56 MeV, p„ = 1.53po, where 
the pressure is already positive and the compressibihty has only one zero (see Fig. 2 b)), is 
identified, as usual, as the critical endpoint, that connects the íirst-order phase transition 
and the crossover regions. At this point we have a second order phase transition once the end 
point of a first-order line is a critical point of the second order. The vàlues of the chemical 
potentials are: /x„ = 304.5 MeV, iJ>d — /J's — 353.3 MeV and = {/lu + + A*s)/3 = 337.0 
MeV. 

In order to get a better insight on the nature of the phase transition in hot and dense 
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FIG. 5: Combined effects of temperature (T) and density (pn/po) in neutron matter in /? - 
equilibrium: a) pressure; b) M„; c) M^. 

matter, we plot in Fig. 5 a) the pressure in the T — p plane. The region of the surface with 
negative curvature corresponds to the region of temperatures and densities where the phase 
transition is first-order. 

We can also illustrate the degree of restoration of chiral symmetry in different flavor 
sectors by plotting the constituent quark masses in the T — p plane. In Fig. 5 b) (see 
also [32]) one can see a clear manifestation of the restoration of chiral symmetry for the 
light quarks that is indicated by the pronounced flattening of the surface with increasing 
temperature and density. Fig. 5 c) shows a more smooth behavior of the strange quark, 
which reflects the well known fact that, as already discussed, chiral symmetry shows a slow 
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tendency to get restored in the strange sector. 



V. MESONS IN COLD QUARK MATTER 

In this section we present our results on the behavior of pseudoscalar mesons in quark 
matter at T = 0, giving special attention to matter in j3 -equilibrium. The presence of 
strange valence quarks is related to a change in different observables. In order to discuss 
the relevance of the strange quark we will show the results for the mesonic behavior also in 
matter without (3 -equilibrium. 

One of our aims is to discuss what can in-medium pseudoscalar meson properties teli us 
about possible restoration of symmetries. If to discuss UaÍX) symmetry it is enough to study 
neutral mesons, this is not sufíicient if one wants to study also chiral symmetry, specially 
in the strange sector. Strangeness enters explicitly in the structure of 77, 77' and of kaons, 
and influences the behavior of pions through the 't Hooft interaction. So, in order to have a 
comprehensive view of restoration of chiral symmetry we have to include the study of pions 
and kaons. 



A. Behavior of pions and kaons 

Hadronic systems created in realistic heavy-ion col·lisions are strongly interacting. Many- 
body excitations may carry the same quantum numbers of the hadrons under investigation 
and, therefore, infiuence hadron properties in medium. Other many-body effects, such as 
the Pauli principle, lead to modifications of hadronic properties, as well. How these effects 
are related with partial restoration of chiral symmetry is a challenging problem in many- 
body physics. As it will be shown, both effects are present in the behavior of the fiavor 
multiplets of pions and kaons. 

Let us analyze the results for the masses of , K~ , K^, and for 7r+, tt*^, 7r~, that are 
plotted as functions of Pn/po in Figs. 6,7 and 9. 

As it was already shown in other works 0,0,01; two kinds of solutions may be found in 
asymmetric matter for pionic and kaonic modes in NJL model, corresponding respectively to 
excitations of the Dirac sea, that are modified by the presence of the medium, and excitations 
of the Fermi sea. Here, in order to appreciate the role of the strangeness degree of freedom. 
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FIG. 6: Pion masses as function of density without (left panel) and with (right panel) (5 - 
equilibrium at T = 0. 

we start by comparing the results obtained in matter with f3 -equilibrium and without, Case 
I. 

In order to get a better insight in the resuhs, we plot the hmits of the Dirac and of the 
Fermi sea continua of qq excitations with the quantum numbers of the mesons under study 
in Figs. 7 and 9 (dashed hnes): cu' = \J + + y^M^j-^-j + is the lower hmit of the 
Dirac continuum, and Uup = \J M"^ + - ^ M'guià) + A^^, uio^ = ^M^T+A^^ - are 
the upper and lower limits of the Fermi sea continuum with Aj the Fermi momentum of 
the quark i. These limits can be obtained by inspection of the expressions for the meson 
propagators (for detalls see 

00)- 

We do not show the corresponding limits for the pions 
because, in the range of densities studied, the pion modes remain outside the continuum. 

Concerning the Dirac sea excitations, we observe the expected splitting between charge 
multiplets: the increase of the masses of K^,K^ and vr" with respect to those of K~,K^ 
and TT"*", respectively, is due to Pauli blocking. At the critical density the antikaons enter 
in the continuum, but, in matter with (3 -equilibrium, they become again bound states 
above po ~ 4po. The difference between the behavior of mesons in matter with and without 
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FIG. 7: Kaon and antikaon masses as function of density at T = 0: Case I (a) and c)) and with 
(3 -equilibrium (b) and d)), -lower limit of the Dirac sea continuum; uJup,í^iow - límits of the 
Fermi sea continuum. 

strange valence quarks (the right and left panels of Figs. 6 and 7, respectively) is more 
evident for kaons than for pions, as expected, since the strange degree of freedom for pions 
only contributes exphcitly through the projector Pab [see ()A3|) ]. and for kaons it contributes 
through its quark structure. The dominant effect is the reduction of the sphtting between 
the kaon and antikaon masses, which is a combination of many-body effects and restoration 
of chiral symmetry (Fig. 7). 

Let's now analyze the other type of solution, which we denoted by the subscript S. Below 
the lower hmit of the Fermi sea continuum of particle-hole excitations there are low bound 
states with quantum numbers of K~ , K'^ ,tt~^ , respectively. These low-energy modes are 
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associated to a first-order phase transition. As a matter of fact, it was shown in 
in the NJL model including vector mesons, where a crossover occurs, the splitting between 
the excitations of the Dirac sea is more pronounced and the low-energy solutions do not 
appear. A general conclusion from our exploration of the high-density region, whether or 
not we have matter in (3 -equilibrium, is that the low-energy modes in this region are poorly 
coUective and, therefore, have little strength, contrary to what happens in the low-density 
region. The pion low-energy mode merges even in Fermi sea continuum just after pn = p'^ 

Since the strangeness degree of freedom is more relevant for kaons than for pions, we will 
analyze the results for kaons in more detail, in particular discussing what happens in Cases 
II and III. Complementary Information concerning the results shown in Figs. 6-7 is given by 
the quark-meson coupling constants displayed in Fig. 8. The coupling constant for kaons 
decreases with density, which is consistent with the increasing of its mass; for antikaons, 
it remains constant when the modes are in the continuum, but, if they get outside of the 
continuum it increases again, which is consistent with the interpretation given above that 
they are bound states in that region of densities. Concerning the low-energy antikaon, it 
becomes less bound as the density increases. 

In Fig. 9 the results for the masses and meson-quark coupling constants are shown for 
Cases II and III. In Case III (matter completely flavor symmetric) all the charged and neu- 
tral kaons and antikaons are degenerated, as expected; its mass increases and its coupling 
constant decreases accordingly. The low-energy antikaon does not exist, which is in agree- 
ment with its origin: a particle-hole excitation due to the flavor asymmetry of matter. The 
behavior for Case II is qualitatively similar to matter in (3 -equilibrium, which is natural, 
since the strange valence quarks only appear at some density. The main differences are that 
the neutral kaons are degenerated with the charged ones, due to isospin symmetry, and the 
upper antikaon gets out the continuum at a higher density than in the case of f3 -equilibrium 
and is less bound, which is due to later appearance and smaller fraction of strange valence 
quarks, in this case. 

Finally, some remarks are in order concerning our results. Below the critical density the 
system is in a mixed phase, therefore, there is no clear deflnition of hadron masses, that 
is why we plotted the masses as dotted lines in this region (they should be understood as 
average vàlues, meaningful only from a qualitative point of view). Moreover, it is clear that 
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FIG. 8: Density dependence of meson-quark coupling constants for T = 0: Case I a) and c); with 
P -equilibrium b) and c). 

a Fermi sea of quarks is certainly not a good description of nuclear matter in the confined 
phase, and the results in that region within NJL models sliould be taken with caution. The 
results for the high-density region are more reliable, since at such densities the quarks are 
supposed to be deconfined and a Fermi sea of quarks is a reasonable description of matter 
in that region. 

A question that can be naturally raised is if our model allows for the possibility of 
kaon condensation. The idea of kaon condensation and the recognition of its relevance 
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42]. In our present approach, 



for astrophysics have been explored since the 80's 
however, the criterion for the occurrence of kaon condensation is not satisfied since the 
antikaon masses are always larger than the difference between the chemical potential of 
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FIG. 9: Kaon and antikaon masses (left panel) and meson-quark coupling constants (right panel) 
as function of density for Case II and Case III at T = 0. In Case II K^, and R'g are degenerated 
with , K~ and Kg respectively. In Case III all kaons are degenerated (dash-dot-dot line). 

strange and nonstrange quarks. Our results should be understood as a starting point for 
further refinements, since it is known that in matter of high-density and low temperatura the 
system might undergo one or more phase transitions. For instance, we have not taken into 
account the effect of color superconductivity that is expected to provide an extra binding 
mechanism. This mechanism is more important for strange quarks, leading to the CFL 
phase 0,0], with kaon condensates, and is supposed to have interesting consequences for 
the structure of neutron stars. This is out of the scope of the present paper, but work in 
this direction is in progress. 

B. Behavior of r/andr/' 

Now, we analyze the results for the masses of rj and rj' . As we can see, the rj' -meson lies 
above the quark-antiquark threshold for p„ < 2.5po and it is a resonant state. After that 
density, the r]' becomes a bound state. This is due to the fact that the limits of the Dirac 
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sea increase with density {uju,d = 2/i„,d and Ug = \JM^ + at p 7^ 0, instead of Ui = 2Mj 
at p = 0). 

Interest in the study of in-medium properties of f], f]' is in part motivated by the conjecture 
that their behavior could give indications of possible restoration of the Ua{^) symmetry, in 
particular their masses could eventually be degenerated. We will show that within our model 
this does not happen, the behavior of rj and rj' essentially reflects the tendency to restoration 
of chiral symmetry in different sectors. 

It has been argued that, in principle, there is no reason why the parameters used to fix 
the model in vacuum should not depend on density or temperature. In a dependence 
on temperature of the coupling constant {gniT) = 5f£)(0)exp[— (T/Tq)^]) was used in order 
to investigate the restoration with temperature of the í/a(1) symmetry, that is explicitly 
broken in this model by the 't Hooft interaction. In the present calculation our parameters, 
including g£, are kept constant, so it would be expected that no indications of restoration of 
Ua{Í) symmetry would be found. 

The masses, plotted in Fig. 10, exhibit a tendency to level crossing but, after the critical 
density, they split again, the splitting being more pronounced in the case of matter without 
strange quarks. This is due to the absence of strange valence quarks in this case and is 
related with the in-medium behavior of the mixing angle, 6, between rj and rj'. As it was 
shown in 3^, above p ~ 3.5po the angle becomes positive and increases rapidly; a similar 
behavior with temperature was found in the framework of the a model Therefore, the 
strange quark content of the mesons changes: at low-density, the rj' is more strange than 
the 1] but the opposite occurs at high-density. So, in a medium without strange quarks the 
f] mass should stay constant in the region where its content is dominated by the strange 
quark. 

Finally, we add some remarks on the approximations made in the study of rj ,1]' in- 
medium. Let us focus again on the projector Pab and Uab ©■ The nondiagonal elements 
that describe the mixing of 7r° — 77 and of 7r° — rj' are proportional to < ququ > — < ÇdQd > 
for Pab and to Juu{Po) — Jdd{Po) for Uab- In the cases here considered of matter with isospin 
asymmetry (matter in f3 -equilibrium and Case I) these quantities are nonzero. Here we 
did the approximation of neglecting them and we checked, by means of a simple estimation, 
that this is a reasonable approximation. 
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FIG. 10: r] and rj' masses as function of density without (left panel) and with (right panel) /3 
-equilibrium at T = 0. tOu ,ijJd - límits of the Dirac continua. 

VI. MESONS IN HOT AND DENSE MATTER 

As it has been discussed in Section IV, the phase transition in hot and dense matter is 
of first-order below ~ 56 MeV. Above we have a crossover. We define the critical 
point as the temperature and density where the pion becomes unbound (Mott transition 
point). Since the model has no confinement, the system is unstable against expansion, 
but, in spite of these drawbacks, we think it is illustrative to plot the meson masses as 
functions of temperature and density. We consider here only the case of neutron matter 
in P -equihbrium. As aheady discussed in Section IV A[ in order to investigate possible 
restoration of symmetries, it is important to study the mesons that have phenomena of 
symmetry breaking at their origin, whether it is chiral or Ua{1) symmetry. As it will be 
shown, the dominant effect found in our results is the restoration of chiral symmetry. 

A first conclusion is that there are differences in the mesonic behavior as compared to the 
zero temperature case discussed in last section, where we have seen that some mesons are 
still bound states in the chiral restored phase. Here, similarly to finite temperature and zero 
density case, the qq threshold for the different mesons is at the sum of the constituent quark 
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FIG. 11: Pion masses as functions of temperature and density. 

masses, so the mesons dissociate at densities and temperatures close to the critical ones. 
A second feature to be noticed is that the mesons that are remnants of Goldstone bosons 
show more clearly the difference between the chiral symmetric and asymmetric phases. The 
shght differences of behavior inside each flavor multiplet are due to many-body effects. It 
can be seen in the diagrams that there is a "hne" separating the regions of the surface with 
different curvatures. This is very clear for the case of pions (Fig. 11) and kaons (Fig. 12), 
and, as shown in js^, this is also apparent for the t]. We may call this line the "Mott circle", 
since it separates the region where the meson are bound states from the region where they 
are in the continuum. 

We notice that in the context of our model, Mott transition is certainly related with 
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FIG. 12: Kaon masses as functions of temperature and density. 

the chiral transition but can not be seen as a mechanism for quark deconfinement. In fact, 
the NJL model is suitable to describe the chiral phase transition but not the deconfinement 
phase transition, since it has no confining mechanism. However, in spite of this drawback, the 
model gives a reasonable description of a quark phase, at high-densities and temperatures, 
wherc intcracting quarks with small constituent mass are supposed to exist. In this phase 
the mesons are unstable resonances, decaying in pairs of quark-antiquark states. 

The behavior of the rj is more involved since this meson is not a Goldstone boson associ- 
ated with spontaneous chiral symmetry breaking, and, in the present model is described as a 
qq resonance, even at zero temperature and density (see Fig. 10). Only at zero temperature 
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and for p„ > S.Opo it becomes a bound state (for more details see |38|). 

Finally, we notice that, as soon as the system heats, and since very low temperatures, 
the low-energy modes stay inside the low continuum and, unlike the zero temperature case, 
they do not become bound states. 



VII. SUMMARY AND CONCLUSIONS 



To conclude and summarize, in the present paper we have investigated phase transitions 
in hot and dense matter and the in-medium behavior of the pseudoscalar mesons, in the 
framework of the SU (3) Nambu-Jona-Lasinio model, including the 't Hooft interaction, 
which breaks the Ua{1) symmetry. 

Three scenarios were considered: i) zero chemical potential and finite temperature, ii) 
zero temperature and finite chemical potential in four types of quark matter and, finally, iii) 
finite temperature and density. Although in all the cases we found (partial) restoration of 
chiral symmetry, different features occur in several observables. 

Concerning case i) we mainly reproduce results obtained by other authors [3^, in order 
to allow for comparison with situations ii) and iii) that represent our original contribution. 
The main feature is the dissociation of mesons at the Mott transition point that occurs when 
the meson masses equal the sum of the masses of their constituents. After that point, which 
sets the corresponding threshold of the qq continuum, the mesons cease to be bound states 
and become resonances. 

New features occur at zero temperature and finite density, in fiavor asymmetric medium 
(ii)). In order to discuss the nature of the phase transition, we plot the pressure and energy 
per particle versus baryonic density. For a suitable choice of the parameters we have a 
mixed phase, which may be interpreted as the system having a hadronic phase with partially 
restored chiral symmetry embedded in a nontrivial vacuum. After the critical density we 
have a quark phase. We notice that in fiavor asymmetric matter the energy of the stable 
hadronic phase, at pcr, is in a region where strange valence quarks are still absent, so SQM 
only could exist in a metastable state. Only for the case of equal number of quarks u ,d ,s we 
found stable SQM, but with a higher energy per particle then atòmic nuclei. Concerning the 
masses of the mesons, there is a splitting between the fiavor multiplets in fiavor asymmetric 
matter and, for kaons and pions, low-energy modes with quantum numbers of vr^, K~ and 
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appear. Our results for kaons in flavor asymmetric matter show that the splitting between 
íaons and antikaons increases with density, which is compatible with experimental results 
2^ that indicate a reduction (enhancement) of kaon (antikaon) production in medium. In 
the high-density region the splitting has a different behavior according to whether there 
are strange quarks present or not. However, although the splitting is reduced in matter 
with strangeness (the smaller splitting is in matter in (3 -equilibrium where there is a larger 
fraction of strange quarks), the upper antikaon becomes more bound in this case, so it 
is likely that the splitting can also be observed. The lower antikaons are more bound in 
the low-density region; these modes do not exist in flavor symmetric matter. Although we 
flnd kaons as bound states in the high-density region, we do not flnd kaon condensation. 
Concerning rj and 77', their masses come closer up to the critical density, but after that point 
they split again, the splitting being more pronounced in the case of neutron matter without 
(5 -equilibrium, which is related with the change of the strangeness content of the mesons. 

In situation iii), hot and dense matter, the phase transition becomes a crossover above the 
critical "end point" T = 56 MeV, p = 1.53po; the system having a mixed phase before that 
point. In this case, the qq continuum turns again to be at the sum of the constituent quark 
masses. By plotting the meson masses we see that there is a line along which the curvature 
of the surfaces changes, indicating the partial restoration of chiral symmetry. Beyond that 
line, the mesons become resonances. 
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APPENDIX A 

In this appendix we give some detalls about the calculation of the effective action and of 
the integrals appearing in the meson propagators, in the vacuum and at flnite temperature 
and density. 
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1. Calculation of the efFective action 



The model Lagrangian (^J can be put in a form suitable for the usual bosonization 
procedure after an adequate treatment of the last term that contains a six quark interaction. 
To obtain a four quark interaction from the six quark interaction we make a shift (çA"g) — > 
{q\°'q)+ < qX"'q >, where < gA"g > is the vacuum expectation value, and contract one 
bilinear {qX°'q). Then the following effective quark Lagrangian is obtained: 

Ceff = qiiYd^ - m)q 

+ Sat[{qX''q){qX'q)]+ Pat[{q^l,X''q){q^l,X'q)], (Al) 

were the projectors Sab , Pab are of the form: 

Sab = gsSab + goDabc < qX^q >, (A2) 
Pab = gsàab - goDabc < qX^q > . (A3) 

The constants Dabc coincide with the SU{3) structure constants datc for a,b,c = (1,2,..., 8) 
and Doab = ~'^^ab, Dqoo = \/f • The hadronization procedure can be done by the integration 
over the quark fields in the functional integral with ()Alj) . leading to the effective action 

2. Integrals 

The polarization operator in Eq. (5) in the momentum space has the form, 

UabiP) = iN, í -^Md \s,{p){X%,{i^,)Sj{p + P)(A''),,(z75)l , (A4) 



. (27r)4- 

where tr^ is the trace over Dirac matrices. 

In the expressions of II*-' in Eq. (5), at T = and p = 0, we have used the following 
integrals: 



d'^p 



(27r)4 {^p2-Mf){{p + Poy-Mf 
Nc p'^d-p Ei + Ej 



471^ Jo EiE, PÍ-{Ei + E,f 
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(A6) 



where Eij = y + is the quark energy. To regularize the integrals we introduce the 
3-dimensional cutoff parameter A. When Pq > Mi + Mj it is necessary to take into account 
the imaginary part of the second integral. It may be found, with help of the ie -prescription 
Pq Pq - that 

- 4vr2 ' Jo E,E, Pi - (E, + E.f + 'l67r (E* + E*) ^^^^ 



with the momentum: p* = ^ {P^ - {Mi - Mjf){Pi - {Mi + Mj)2)/2Po and the energy: 
El, = vW+mJ. 

Concerning the calculation of the propagators of the diagonal mesons 7r°, 77 and rj', the 
projector Pab and the polarization operator Hat (see Eq. (|7j)) are matrices that, in the case 
< QuQu >=< (ídÇd >, have the nonvanishing elements: 



^33 = gs + go <qs Qs >, (A8) 
2 

^00 = gs- -^go {< ququ > + < qdqd > + < qs qs >) , (A9) 

^88 = gs + ^go {2 < q^qu > +2 < qdqd > - < qs qs >) , (AlO) 

P08 = Pso = ^j^9d {< ququ> + < qd qd> -2 < qs qs >) , (All) 



and 



noo(Po) = ^ [Juu{Po) + Jdd{Po) + Jss{Po)] , (A12) 

n88(Po) = \ [Juu{Po) + Jdd{PQ) + 4J,,(Po)] , (A13) 

no8(Po) = n8o(Po) = ^ [Juu{Po) + Jdd{Po) - 2Jss{Po)] , (AM) 



where 



UPo) =^{l\- ^mPoij . (A15) 

At finite temperature and with two chemical potentials, /i, and /Xj, the integrals I[ , I2 , 12 
take the form: 
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I'2ÍPo,T,IJ,i,IJ,j) = -iVc 



d'p 


1 


1 


(27r)3 


2Ei {Ei + Po - 
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1 




2Ei {E, -Po + (/i, - 
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' 2E^ {Ej - Po + (/ií 


-f^jW-Ef 


1 


1 





Í2 



n- 



11 . 

2Ej{E,+Po-{i2,-i2,)y-E! ' 

where nf are the Fermi distribution functions ()14p. defined in Section III. 
For the case i = j, with imaginary part, we have the expression 

Nc p f p'rfp 1 
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